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ABSTRACT 

We investigate the dissociation equilibrium of H 2 in very cool, helium-rich white dwarf at¬ 
mospheres. We present the solution of the non-ideal chemical equilibrium for the dissociation of 
molecular hydrogen in a medium of dense helium. We find that at the photosphere of cool white 
dwarfs of T^s = 4000 K, the non-ideality results in an increase of the mole fraction of molecular 
hydrogen by up to a factor of ~ 10, compared to the equilibrium value for the ideal gas. This 
increases the H 2 — He CIA opacity by an order of magnitude and will affect the determination of 
the abundance of hydrogen in very cool, helium-rich white dwarfs. 

Subject headings: dense matter - stars - atmospheres - stars: white dwarfs 


1. Introduction 

Several very cool white dwarfs with suspected 
Tes ^ 4500 K have been discovered recently (Far- 
ihi 2005; Kilic et al. 2005; Gates et al. 2004; Op- 
penheimer et al. 2001; Harris et al. 2001; Hodgkin 
et al. 2000; Ibata et al. 2000; Harris et al. 1999). 
Most of them are thought to posses helium-rich 
atmospheres with an very high He/H ^ 10^ ratio 
(Bergeron et al. 2005; Kilic et al. 2005; Gates et 
al. 2004; Bergeron & Leggett 2002; Bergeron 2001; 
Oppenheimer et al. 2001; Hodgkin et al. 2000). In 
most cases, however, current atmosphere models 
fail to reproduce the observed spectra and pho¬ 
tometry of these peculiar stars. The reason, and 
there may be more than one, for this shortcoming 
of the models is presently unknown. However, cur¬ 
rent models predict extreme physical atmospheric 
conditions for such stars, reaching densities of up 
to 2 — 3 g/cm^. Under these conditions, the mostly 
ideal gas constitutive physics used in published at¬ 
mosphere models is demonstrably inadequate. A 
careful look at the dense matter effects on the 
equation of state, chemistry, opacities, and radia¬ 
tive transfer is necessary to compute physically re¬ 
alistic models of these stars. Several of these ef¬ 


fects have been studied previously, such as refrac¬ 
tive radiative transfer (Kowalski & Saumon 2004), 
the effects of fluid correlations on He“ free-free 
and He Rayleigh scattering (Kowalski et al. 2005; 
Iglesias et al. 2002), and the ionization of warm, 
dense helium (Kowalski et al. 2005; Bergeron et 
al. 1995). In this contribution we present an addi¬ 
tional correction that arises in the dense fluid: The 
solution for the dissociation of molecular hydro¬ 
gen in dense fluid helium, in the limit He/H >> 1. 
The relative importance of these corrections varies 
considerably, even more so when they are com¬ 
bined. As several more dense matter effects re¬ 
main unexplored, it is premature to ponder their 
implications for the analysis of the coolest white 
dwarfs known and whether they will result in mod¬ 
els that reproduce the data. Nonetheless, incor¬ 
porating adequate constitutive physics in atmo¬ 
sphere models is a necessary step to reach a proper 
understanding of these peculiar stars. 

We introduce non-ideal effects into the equilib¬ 
rium dissociation of molecular hydrogen through a 
modification of the chemical potentials of H i and 
H 2 (section 2). We find that the strong interac¬ 
tions in the dense, helium-rich atmosphere results 
in a significant decrease in the dissociation frac- 
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tion of molecular hydrogen, with a corresponding 
change in the H 2 — He Collision-Induced Absorp¬ 
tion (CIA) opacity, which is linear function of uhs- 
In section 3, we illustrate the impact of the inter¬ 
actions on the H 2 /H i ratio on a sequence of white 
dwarf atmosphere models with Teff = 4000 K, a 
gravity of log g = 8 (cgs), and a homogeneous 
composition of He/H = 10^, 10^, and 10®, where 
He/H is the number abundance ratio. 


2. The dissociation equilibrium of molec¬ 
ular hydrogen in a dense fluid 

2.1. Theoretical approach 

The condition for chemical equilibrium (at a 
given density and temperature) for the dissocia¬ 
tion reaction: 

H2v^2Hi (1) 

is given by (Cox & Giuli 1968) 

~ 2mhi = 0, (2) 


where /ii is a chemical potential of the species i 
expressed as 


Mi = Eo^i -I- ksTln 


Zi(27rmifcBr)®/2 


Mi 


(3) 

In the above equation, is the Boltzmann con¬ 
stant, h the Planck constant, T the temperature, 
Aoy the ground state energy, m the number den¬ 
sity, Zi the unperturbed internal partition func¬ 
tion, and mi is the mass. The first two terms on 
the r.h.s. of equation (3) represent the ideal con¬ 
tributions of translational and internal degrees of 
freedom and is the non-ideal contribution 

to the chemical potential arising from the interpar¬ 
ticle interactions in the fluid. Setting = Q 

we recover the standard Saha equation for disso¬ 
ciation of molecular hydrogen: 


r,id ^ ^ 

"-Hi ^Hi 


mu, 




[2TTm‘^^kBT \ 


3/2 


^DojkBT 


(4) 


If we define the quantity A/ as 

A/ = - SmST”*"'’ (5) 

the non-ideal equilibrium equation can be written 
in the following form 

1,2 i3/2 

mn^n- ^(Do-Ai)/kBT 

2TTrn^-JiBT 

= ( 6 ) 

Comparing (4) with (6), we see that the non¬ 
ideal effects on the dissociation equilibrium can be 
interpreted as a change in the dissociation energy 
by a value of A/. For simplicity, we will follow 
this interpretation hereafter. We emphasize that 
this description of the non-ideal contribution to 
the chemical equilibrium (Eq. 6) is identical to the 
occupation probability formalism of Hummer & 
Mihalas (1988) if there is only one bound state in 
the partition function (e.g. low temperature Hi) 
or if the is the same for all bound states 

j of the species i = df{V,T,nij)/dnij 

remains constant in Eq. (2.17) of Hummer & Mi¬ 
halas (1988)). In both cases, the occupation prob¬ 
ability can be factored out of the partition function 
and written as 0 ^^ 

In the atmosphere of cool white dwarfs, hydro¬ 
gen exist mostly as H 2 and Hi, and the H 1 /H 2 
ratio is governed by reaction (1) only. For a given 
temperature, density p, and the atmosphere com¬ 
position y = He/H, the number densities of H 2 
and H i are: 


_ ^ 


riBi = 


~1 + VI + 8/3ntot 
Ad ^ 

7 IH 2 I ’ 


where 


ntot = + 2nH2 = 


P 


TOh + yWHe 
refers to the hydrogen species only. 


(7) 

( 8 ) 

(9) 


where Dq = 2i?o,Hi—ifo,H 2 = 4.478 eV is the disso¬ 
ciation energy of the isolated hydrogen molecule. 
Even for trace species, like Hi or H 2 in dense he¬ 
lium, the which arise from interactions 

with the atoms are not negligible and in principle 
should be comparable in magnitude to 


2.2. Computation of the H 2 dissociation 
equilibrium 

The non-ideal contributions to the chemical po¬ 
tentials of HI and H 2 were obtained through the 
numerical solution of the Ornstein-Zernike equa¬ 
tion in the Percus-Yevick (PY) approximation 


2 




r (a.u.) 

Fig. 1.— The H — He and H 2 — He pair interaction potentials. The dashed line represents the radius 
(horizontal axis) of a sphere whose volume equals the effective volume of two H atoms colliding with a 
kinetic energy of ffesT (see text). 


(Martynov 1992). For the H — He interaction we 
use the pair potential of Shalabi et al. (1998), and 
for the H 2 — He interaction the pair potential given 
by Ree (1983) (Figure 1). Both potentials are 
from ab initio quantum mechanical calculations 
and are in good agreement with the independent 
calculations of Tang & Yang (1990) for H — He 
and of Tao (1993) and Shafer & Gordon (1976) 
for H 2 — He. As we consider a helium-dominated 
mixture (He/H ^ 10^), the H — H, H 2 — H, and 
H 2 — H 2 interactions can be neglected. 

High pressure experiments have shown that ab 
initio pair potentials are too repulsive to describe 
dense systems where N-body effects become im¬ 
portant (Nellis et al. 1984; Ross et al. 1983). The 
softening of the pair potentials at high densities 
can only be quantified experimentally or, alterna¬ 
tively, estimated with N-body quantum mechan¬ 
ical calculations. Since neither are available for 
mixtures of trace hydrogen in helium, we resort to 
ab initio potentials. The net effect on the dissocia¬ 
tion equilibrium depends on the relative softening 
of the potentials (Eq. 5) and is therefore less sen¬ 
sitive to N-body effects than the individual poten¬ 
tials. We also calculated the chemical potentials in 
the Hypernetted Chain approximation and found 


them to agree within 5% with the PY values up 
to 2g/cm^. Since the PY approximation is better 
suited for short range potentials such as the ones 
we use here, we estimate that our PY calculations 
are reliable up to at least 2g/cm^. 

For the internal partition functions Zi, we use 
expressions for the electronic ground state of the 
unperturbed hydrogen molecule accounting for 
the vibrational/rotational excitations (Huber & 
Herzberg 1979), and set Zhi = 2 for hydrogen 
atom. This approximation is justified as the elec¬ 
tronic excitation energies of both species are large 
and for temperatures of a few thousands degrees 
the populations of the electronic excited levels 
are extremely small. However, there is significant 
thermal excitation of the rotational and vibra¬ 
tional levels of H 2 and the effect of the dense 
medium on must be considered. Since the 
molecule does not have spherical symmetry, it is 
possible that its rotation modes will be hindered 
by interactions with neighbors at very high den¬ 
sity. Furthermore, the energies of the vibrational 
levels, which are spaced by a mere 0.54 eV, could 
also be significantly shifted by these interactions. 
Both of these effects are discussed by Saumon & 
Chabrier (1991) and can be neglected under the 
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Fig. 2.— (a) Variation of the effective dissociation energy Dq — A/ as a function of the density of helium 
for temperatures of 3000 K (solid line), 4000 K (dotted line), and 5000 K (dashed line), (b) The density p at 
the photosphere of cool white dwarf atmosphere models of = 4000 K, log 5 = 8 (cgs), and various He/H 
compositions. 


present conditions. The use of the unperturbed 
rotational and vibrational levels of H 2 is thus jus¬ 
tified in the calculation of Zhs ■ 

[hb] 

More importantly, the excited rotational and 
vibrational levels of the H 2 molecule may be differ¬ 
entially affected by the interactions. This relates 
to the general problem of the cutoff of the sum over 
states in the internal partition function and is the 
subject of an extensive literature. This truncation 
of the partition function sum is often described by 
an occupation probability formalism (Hummer & 
Mihalas 1988) and in the case of interactions be¬ 
tween neutral particles, is generally described in 
the context of the volume occupied by each bound 
state. In the case of a diatomic molecule, the vol¬ 
ume occupied increases with the rotational quan¬ 
tum number (centrifugal stretching) and with the 
vibrational quantum number (larger amplitude of 
vibration). The first effect is very small and can 
be safely neglected. We estimate the vibrational 
stretching by considering the change in the equilib¬ 
rium internuclear separation of the molecule with 
the vibrational quantum number. Because of the 
anharmonicity of the potential, higher vibrational 
levels correspond to larger equilibrium separations 
and a larger average molecular volume. Based on 
the potential curve of the H 2 molecule (Kolos & 


Wolniewicz 1965) we find that for T < 6000 K and 
densities of up to 2 g/cm^, the vibrational excita¬ 
tion is largely limited to the lowest five levels and 
the partition function is reduced by a few percent 
when the larger volume of the excited states is 
taken into account. This is a small effect which 
can be neglected in view of the other uncertainties 
in the model such as the N-body effects on the 
interaction potentials. Thus, the use of the inter¬ 
nal partition function of the isolated H 2 molecule, 
with all states being affected identically by the in¬ 
teractions with He (Eq. 6 ) is a very good approx¬ 
imation. 

3. Results and Discussion 

3.1. The H 2 dissociation equilibrium in 
fluid helium 

In general, the A/ defined by (5) is a function 
of the density of helium and the temperature (Fig. 
2a), but is independent of He/H ratio if H is a 
trace species (He/H ^ 10^). For the conditions 
in white dwarf atmospheres, the temperature de¬ 
pendence is weak. The change in the dissociation 
energy A/ is negative, making H 2 more stable in 
dense helium. This may be qualitatively under¬ 
stood by comparing the effective volume occupied 
by one H 2 molecule to that of two H atoms in He. 
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Table 1 

The change in the dissociation energy AI (eV) for representative densities and 

TEMPERATURES. 


P (g/cm^) 

r(K) : 2000 

4000 

6000 

8000 

0.001 

-0.001 

-0.001 

-0.001 

-0.001 

0.010 

-0.007 

-0.010 

-0.011 

-0.012 

0.100 

-0.083 

-0.108 

-0.120 

-0.127 

0.250 

-0.242 

-0.297 

-0.323 

-0.334 

0.500 

-0.579 

-0.668 

-0.697 

-0.702 

0.750 

-1.006 

-1.075 

-1.088 

-1.070 

1.000 

-1.471 

-1.494 

-1.469 

-1.414 

1.500 

-2.011 

-2.094 

-1.963 

-1.825 

2.000 

-2.189 

-2.387 

-2.192 

-1.988 


Lets assume that these effective volumes are the 
spaces around each of the particles where the en¬ 
ergy of their interaction with He atoms is greater 
than the average thermal kinetic energy of the par¬ 
ticles IfcsT. This radius is the classical distance 
of closest approach in a collision. On Figure 1, 
we plotted the radius of a sphere whose volume 
equals the effective volume of two hydrogen atoms, 
as a function of ksT (dashed curve). The corre¬ 
sponding radius for the hydrogen molecule is rep¬ 
resented by the H2 — He potential (solid curve). 
The effective volume of H2 is smaller than that 
of two H atoms for fc^T ^ leV. Since the exclu¬ 
sion of a greater volume results in a decrease of 
the entropy of He, H2 is more stable in dense he¬ 
lium than two H atoms. In dense hydrogen, we 
have the opposite situation where H2 is less stable 
at high density (leading to pressure dissociation 
of H2) because H2-H2 is more repulsive than the 
sum of the H-H and H-H2 interactions (Saumon 
& Chabrier 1991). 

Numerical values for AI obtained from our 
model are given in Table 1. Together with rela¬ 
tions (6)-(9) this table allows for an easy imple¬ 
mentation of the non-ideal abundances of atomic 
and molecular hydrogen into white dwarf atmo¬ 
sphere codes. 

The non-ideal recombination of H becomes sig¬ 
nificant when —AI/ksT ^ 1. For an atmosphere 
with Teff = 4000 K, this will occur at the photo¬ 
sphere if p ^ 0.3g/cm^, a density easily achieved 
in He-rich models. As Dg — AI increases with in¬ 


creasing helium density, a lower dissociation frac¬ 
tion will result. A strong deviation from the ideal- 
gas abundance of H2 is expected. 

3.2. The H 2 dissociation equilibrium in 
He-rich white dwarf atmospheres 

We have applied the non-ideal correction to 
the H2 dissociation equilibrium to very cool white 
dwarf atmospheres. To show the effect, we solved 
for the chemistry and opacity along fixed T — p 
white dwarf atmosphere profiles of T^s = 4000 K, 
a gravity of log 5 = 8 (cgs), and various He/H ra¬ 
tios. Figure 3 shows the atmospheric density pro¬ 
files. These profiles were obtained with a code that 
solves for static, plane-parallel, and LTE white 
dwarf atmospheres in thermal and hydrostatic 
equilibrium, accounting for refraction (Kowalski 
& Saumon 2004). The resulting new abundances 
of Hi and H2 in our models are given on Figure 4 
(left panels). As expected, the amount of molec¬ 
ular hydrogen increases significantly throughout 
the atmosphere for He/H ^ 10^, and H2 can be¬ 
come the dominant hydrogen species at the pho¬ 
tosphere. 

The most important sources of opacity in 
helium-rich white dwarf atmospheres are He“ 
free-free, Rayleigh scattering, and H2 — He CIA 
(Hansen 2004). Several physical and chemical ef¬ 
fects that alter these opacity sources have already 
been discussed: the change in the free-free opacity 
and the Rayleigh scattering of helium, as caused 
by the strong correlations in the dense fluid (Igle- 
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Fig. 3.— The density profiles of white dwarf atmosphere models of T^s = 4000 K, log g = 8 (cgs), and He/H 
composition of He/H = 10®, 10^ and 10^ (from top to bottom, respectively). 


sias et al. 2002), the change in the number density 
of free electrons (Kowalski et al. 2005; Bergeron et 
al. 1995), the presence of heavy elements (Berg¬ 
eron 2001), and the formation of trace species like 
HeJ (Malo et al. 1999) and HeH+ (Harris et al. 
2004). The effect on non-ideal chemistry on the 
H2 dissociation on the opacity is shown in the 
right hand panels of Figure 4. The increase in 
number density of molecular hydrogen results in 
an increase of the H2 — He CIA opacity. This 
follows from the linear dependence of H2 — He 
CIA opacity on riHa- On the other hand, the 
H“ bound-free and free-free opacities are reduced, 
due to the decrease in the abundance of atomic 
hydrogen. Based on the atmospheric structures 
used here, the effects of the non-ideal chemistry of 
H2 in dense helium are maximal for He/H ~ 10®. 
This arises from a competition between the need 
for a high He/H ratio to increase the non-ideal ef¬ 
fects by increasing the density at the photosphere 
(Fig. 2b) and the need for a high enough hydrogen 
content in the atmosphere so that H2 (or H) con¬ 
tributes to the total opacity (Fig. 4, right hand 
side). 

The impact of the non-ideal dissociation equi¬ 
librium of H2 on the synthetic spectrum of a He- 


rich white dwarf model is shown in Figure 5. For 
a white dwarf atmosphere model of Tgg = 4000K, 
log g = 8, and He/H = 10®, the rise in abun¬ 
dance of H2 increases significantly the opacity 
in the infrared, causing a redistribution of the 
flux toward shorter wavelengths. This effect on 
the spectrum of cool white dwarfs is largest for 
Teff = 4000K —4500K, where H2 is partially disso¬ 
ciated. At lower effective temperatures hydrogen 
exists mostly in molecular form and the effect of 
the non-ideal dissociation equilibrium on the spec¬ 
trum vanishes at Tgg ~ 3000K. 

The determination of the He/H composition of 
very cool He-rich white dwarfs depends mostly on 
the relative He — H2 CIA, and He“ free-free opac¬ 
ities (Fig. 4, right hand side). Since the non-ideal 
recombination increases the CIA opacity, this ef¬ 
fect, taken by itself, will result in a higher value 
for the He/H ratio inferred from a given stellar 
spectrum. On the other hand, current calcula¬ 
tions for the ionization fraction of dense helium 
and the low-frequency behavior of the He“ free- 
free opacity are uncertain and the He“ free-free 
opacity could be underestimated by as much as 
~ 2 —4 orders of magnitude (Kowalski et al. 2005). 
It is therefore too early to draw conclusions about 
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log Tr a (/xm) 

Fig. 4.— Left panels: The number density of atomic (solid) and molecular (dotted) hydrogen with (thick 
lines) and without (thin lines) the non-ideal chemical equilibrium along the atmosphere profiles of Figure 3. 
Right panels: Corresponding contributions to the photospheric opacity: H2 — He CIA (solid), He“ free-free 
of Iglesias, Rogers, & Saumon (2002) (long dashed), H“ bound-free and free-free (dotted) and Rayleigh 
scattering (short dashed). 
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Fig. 5.— The synthetic spectrum for the white dwarf atmosphere models of Tgg = 4000 K, log g = 8 (cgs), 
and He/H = 10^ with (thick line) and without (thin line) the non-ideal dissociation equilibrium of H2. 


the atmospheric composition of the coolest white 
dwarfs. However, we emphasize that, as the chem¬ 
ical equilibrium (6) is not affected by the number 
of free electrons in the atmosphere, our solution 
for the non-ideal abundance of H2 will not be af¬ 
fected by weak ionization of He, and the results 
are limited only by the validity of the H2 — He and 
H — He interaction potentials. 

The non-ideal dissociation equilibrium in pure 
hydrogen atmospheres is discussed in Saumon & 
Jacobson (1999). As hydrogen is much more 
opaque than helium, the density at the photo¬ 
spheres of these stars is much smaller than for 
helium-rich atmosphere stars of the same effec¬ 
tive temperature. The non-ideal effects are there¬ 
fore weaker and become important only in stars of 
lower effective temperatures and/or higher gravity, 
where high densities are achieved. Saumon & Ja¬ 
cobson (1999) show that the non-ideality matters 
in pure-hydrogen atmospheres with Teff ^ 2500K. 
White dwarfs with hydrogen-rich atmospheres of 
such low effective temperatures have not yet been 
identified. 


4. Conclusions 

Recent discoveries of a number of cool white 
dwarfs with peculiar spectral energy distribution 
represent a challenge in the modeling of very cool 
white dwarf atmospheres, as the spectra of these 
stars cannot be fitted with existing models. We 
believe that the main reason of this shortcoming of 
the models is the poorly explored, extreme phys¬ 
ical regime found inside these atmospheres. We 
presented a correction to the abundance of H2 
that arises from the strong interaction of hydro¬ 
gen molecules and atoms in a dense, fluid helium 
medium. We have found that in the dense, helium- 
rich, cool white dwarf atmospheres the formation 
of H2 is more favorable than in the ideal gas de¬ 
scription. For white dwarfs of = 4000 K and 
He/H ^10^ the abundance of molecular hydrogen 
increases by an order of magnitude, with a corre¬ 
sponding increase in the H2 — He CIA opacity by 
the same factor. This improvement is a new, sig¬ 
nificant effect that must be included in realistic 
modeling of very cool white dwarf atmospheres. 

We have shown that the non-ideal effects af¬ 
fect strongly the abundances of even trace species 
in the atmosphere. This strongly suggest the ne- 




cessity of revising the abundances of other trace 
species with signihcant opacity for similar effects. 
We expect that the study of their abundances and 
absorption processes in fluid helium will signifi¬ 
cantly improve our understanding of the atmo¬ 
spheric physics, composition, and evolution of the 
oldest and coolest white dwarfs. 

I thank D. Saumon for useful discussions and 
the referee, P. Bergeron, for suggestions that im¬ 
proved the clarity of the manuscript. This research 
was supported by the United States Department 
of Energy under contract W-7405-ENG-36. 
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